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ABSTRACT 

We present deep NH3 observations of the L1495-B218 filaments in the Taurus molecular 
cloud covering over a 3 degree angular range using the K-band focal plane array on the 100m 
Green Bank Telescope. The L1495-B218 filaments form an interconnected, nearby, large complex 
extending over 8 pc. We observed NH 3 (1,1) and (2,2) with a spectral resolution of 0.038 km/s 
and a spatial resolution of 31". Most of the ammonia peaks coincide with intensity peaks in dust 
continuum maps at 350 pm and 500 pm. We deduced physical properties by fitting a model 
to the observed spectra. We find gas kinetic temperatures of 8 — 15 K, velocity dispersions of 
0.05 — 0.25 km/s, and NH 3 column densities of 5xl0 12 — lxlO 14 cm -2 . The CSAR algorithm, 
which is a hybrid of seeded-watershed and binary dendrogram algorithms, identifies a total of 
55 NH 3 structures including 39 leaves and 16 branches. The masses of the NH 3 sources range 
from 0.05 M 0 to 9.5 M 0 . The masses of NH 3 leaves are mostly smaller than their corresponding 
virial mass estimated from their internal and gravitational energies, which suggests these leaves 
are gravitationally unbound structures. 9 out of 39 NH 3 leaves are gravitationally bound and 
7 out of 9 gravitationally bound NH 3 leaves are associated with star formation. We also found 
that 12 out of 30 gravitationally unbound leaves are pressure-confined. Our data suggest that 
a dense core may form as a pressure-confined structure, evolve to a gravitationally bound core, 
and undergo collapse to form a protostar. 

Subject headings: ISM: clouds — ISM: molecules — radio lines: ISM — stars: formation 

1. INTRODUCTION 
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Filaments are common structures in molecular 
clouds (e.g. Ridge et al. 2006; Andre et al. 2010, 



2013; Men’shchikov et al. 2010). They are im¬ 
portant structures in star formation because most 
dense cores form within filaments and stars form 
within dense cores. Studies of isolated dense cores 
suggest that stars form mainly as a result of grav¬ 
ity within dense cores, but the formation of fila¬ 
ments and the evolution from filaments to dense 
cores are still unclear. Recent theoretical stud¬ 
ies of molecular clouds suggest that a filament 
may form by colliding flows of turbulent media 
and some dense parts further undergo contrac¬ 
tion by either self-gravity or ram pressure of col¬ 
liding flows to make dense cores (e.g. Balsara 
et al. 2001; Padoan et al. 2001; Klessen et al. 
2005; Vazquez-Semadeni et al. 2005; Gomez et 
al. 2007; Nakamura & Li 2008; Gong & Ostriker 
2009, 2011; Heitsch et al. 2008, 2009; Chen et al. 
2014). Since filaments are formed through collid¬ 
ing flows in this scenario, young and less massive 
dense cores within the filaments are gravitation¬ 
ally unbound, while evolved and massive dense 
cores tend to be gravitationally bound. Thus, ob- 
servationally studying the physical state of dense 
cores within filaments with respect to their evo¬ 
lutionary stage is an essential step to understand 
how low-mass star formation proceeds within fila¬ 
ments. 

The Taurus molecular cloud is one of the best 
testbeds for studies of both filaments and dense 
cores because it is nearby (D ~ 140 pc, Loinard 
et al. 2007, Torres et al. 2009) and forming low- 
mass stars (no O and B stars) in a relatively iso¬ 
lated mode (with little spatially confusion, see re¬ 
view by Kenyon, Gomez, & Whitney 2008). The 
L1495-B218 filaments are prominent in the Tau¬ 
rus molecular cloud extending over 3 degrees on 
the sky. They are known to contain cores that 
range from chemically young starless cores (e.g. 
L1521B and L1521E; Hirota et al. 2004, Tafall & 
Santiago 2004) to more evolved protostellar cores 
(in B7/L1495; see Rebull et al. 2010). In addi¬ 
tion, this region is well mapped in dust extinc¬ 
tion (Schmalzl et al. 2010), dust continuum emis¬ 
sion (Herchel Space Observatory , Palmeirim et al. 
2013), dust polarization (Chapman et al. 2011), 
low density molecular gas tracers ( 12 CO and 13 CO 
Goldsmith et al. 2008; C 18 0 Onishi et al. 1996; 
Hacar & Tafalla 2013.), and a dense gas tracer 
(N 2 H+, Hacar et al. 2013). These maps provide 
basic information on the structure of the filaments 


and the positions of some of the densest cores. 

Hacar & Tafalla (2013) observed the L1495- 
B218 filaments in C 18 0 J = 1—0 and N 2 H + J = 
1—0. C 18 0 J = 1—0 is a moderate density (n e ff < 
10 3 cm -3 at 10 K^gas tracer and is a good tracer 
with which to study velocity structure of moder¬ 
ate density regions since its optical depth is much 
smaller than those of 12 CO and 13 CO. Using C 18 0, 
they estimated the LSR velocities and found that 
the L1495-B218 filaments are not a single coherent 
velocity structure but a complex group of at least 
35 interwoven filaments. However, C 18 0 does not 
trace dense cores because CO is depleted onto dust 
grains at a high density and low temperature (e.g. 
Kuiper et al. 1996; Willacy et al. 1998; Krae- 
mer et al. 1999; Caselli et al. 1999; Bacmann 
et al. 2002; Bergin et al. 2002; Tafalla et al. 
2002; Pontoppidan 2006 ; Ford & Shirley 2011). 
To find dense cores, the authors used N 2 H + J = 
1 — 0 , which is a dense gas tracer with an effec¬ 
tive excitation density above 10 4 cm ~ 3 (Reiter et 
al. 2011), but the densities that N 2 H + traces are 
much higher than those traced by C 18 0, which 
suggests that N 2 H + does not trace the complete 
population of both young and evolved dense cores 
within the filaments. In order to identify the com¬ 
plete population of dense cores and probe the tran¬ 
sition between filaments and dense cores, obser¬ 
vations that trace gas with densities intermediate 
between those traced by C 18 0 and N 2 H + are re¬ 
quired. 

Ammonia is an excellent probe that has an 
effective excitation density intermediate between 
C 18 0 and N 2 H + . Ammonia starts to form at an 
early stage of prestellar core evolution and be¬ 
comes brightest during the later stage of prestel¬ 
lar core evolution (Suzuki et al. 1992). NH 3 
traces densities >10 3 cm -3 (Flower et al. 2006; Di 
Francesco et al. 2007). Owing to the NH 3 chem¬ 
istry, the ammonia abundance steeply increases as 
a molecular cloud core becomes denser and it may 
be the most abundant molecule after CO that is 
easily observed within dense cores (Tafalla et al. 
2004; Aikawaetal. 2005). Ammonia is also a good 
tracer for probing the physical condition of fila¬ 
ments. The inversion transitions of ammonia are 


1 n e fi is the effective excitation density, which is needed to 
produce = 1 K line (Evans 1999; Shirley 2015). See 
Reiter et al. (2011) for a summary of effective excitation 
densities. 
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known to provide an accurate measure of gas ki¬ 
netic temperature along the line-of-sight (e.g. Ho 

6 Townes 1983; Rosolowsky et al. 2008) since 
a transition between meta-stable energy levels is 
possible only through collisions. If the signal to 
noise ratio is large enough, the uncertainty of a 
kinetic temperature measured using the inversion 
transitions of ammonia may be as low as 0.1 K. 
An accurate measurement of the kinetic temper¬ 
ature also permits an accurate determination of 
non-thermal motions, making ammonia an excel¬ 
lent tracer for investigating kinematics. 

We observed the L1495-B218 filaments in NH 3 
(1,1) and (2,2). The complete maps spans over 
three degrees on the sky and is the largest NH 3 
map published to date. In this paper, we present 
observations of filaments and dense cores traced 
by the ammonia emission. We analyze the physi¬ 
cal properties of nested structure identified in the 
NH 3 emission. The kinematics will be explored 
in a subsequent paper. The layout of this paper 
is as follows: §2 briefly introduces our mapping 
of the L1495-B218 filaments in NH 3 and the data 
reduction procedures. In §3 we present the ba¬ 
sic physical properties of the L1495-B218 filaments 
and the dense structures identified by the CSAR 
algorithm within the filaments. In §4 we discuss 
physical states of the dense structures identified by 
the CSAR algorithm. Finally, in §5 we summarize 
our results. 

2. OBSERVATION & DATA REDUC¬ 
TION 

2.1. Mapping of the L1495-B218 Fila¬ 
ments in NH 3 (1,1) and (2,2) 

We observed the L1495-B218 filaments in the 
Taurus cloud using the 100 meter Robert C. Byrd 
Green Bank Telescope (GBT12A295)[^] Observa¬ 
tions were carried out from April 2012 to Febru¬ 
ary 2013 in 13 shifts totalling 66 hours. We used 
the K-band focal plane array (KFPA), which has 

7 beams each with dual polarization in a hexag¬ 
onal arrangement (for technical properties of the 
KFPA, see https://science.nrao.edu/facilities 
/gbt/observing/GBTog.pdf). We used the 7+1 


2 The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under coopera¬ 
tive agreement by Associated Universities, Inc. 


mode of the KFPA, which allows the center beam 
to observe two frequencies simultaneously while 
the other 6 beams observe only a single frequency. 
The bandwidth of each frequency window is 50 
MHz. For 7 beams, we set the rest frequency to 
be 23.702 GHz, which between the strongest hy- 
perfine component of NH 3 (1,1) (23.6944955 GHz; 
Lovas & Dragoset 2003) and NH 3 (2,2) (23.7226333 
GHz; Lovas & Dragoset 2003). For the central 
beam we selected a second frequency to include 
the CCS Nj = 2i —>■ lo transition (22.344033 GHz; 
Yamamoto et al. 1990); these results will be re¬ 
ported in a separate paper. The spectrometer is 
configured to have 16384 channels across 50 MHz, 
yielding a resolution of 3.05 KHz (0.0386 km s -1 
at 23.705 GHz). Since the typical thermal broad¬ 
ening of starless cores is < 0.2 km s -1 (e.g. Jijina 
et al. 1999; Caselli et al. 2002; Rosolowsky et al. 
2008), the line profiles of NH 3 are well resolved in 
this setup. 

We observed the L1495-B218 filaments in 20 
rectangular regions that are arranged to cover high 
extinction regions (Schmalzl et al. 2010) while 
minimizing the total observation time. Each rect¬ 
angular region is 20' by 6' or 6' by 20' oriented in 
Galactic coordinates. Data for each rectangular 
region are obtained by on-the-fly mapping of 25 
rows 20' long with a row spacing of 15.12". The 
integration time per row is 150 seconds. The re¬ 
sulting integration time for a direction observed 
with all 7 beams is ~30 seconds, the total obser¬ 
vation time for each rectangular region is 1 hour 
and 20 minutes. The reference OFF positions are 
picked to be where the extinction is low enough to 
suggest that they would have negligible NH 3 emis¬ 
sion, and be close to each rectangular region. The 
average of two observations of the OFF position 
made before and after mapping each rectangular 
region is subtracted from each on-source observa¬ 
tion. We also carried out “daisy” pattern obser¬ 
vations to improve the signal-to-noise ratio (SNR) 
toward 29 bright NH 3 regions that we found from 
the mapping of rectangular regions. Each daisy 
pattern observation was made with a radius of 
2.8' and 22 cycles following a rose curve. The 
daisy pattern observation completely covers the 
inner 2.8' radius with all 7 beams for an average 
integration time of 38 seconds. The total inte¬ 
gration time for each daisy pattern is 20 minutes. 
We updated the pointing of the telescope using 
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the quasar 0403+2600 between every mapping of 
a rectangular region or every three daisy patterns. 
The pointing corrections were usually less than 
10", compared to the beam size of the GBT at 
23 GHz of 31". 

Data reduction was performed using the GBT 
KFPA pipeline. The data are regridded into regu¬ 
lar spaced pixels using AIPS and the gain of the 
7 beams are calibrated through observing Venus 
and Jupiter in every observing shift. As a dou¬ 
ble check on the absolute flux calibration, we 
also observed the peak dust continuum position 
of the L1489PPC starless core (a = 04:04:47.6, 
<5 = +26:19:17.9, J 2000.0) every observing shift 
(Young et al. 2004; Ford & Shirley 2011). The 
calibrations using Venus and Jupiter and the cal¬ 
ibration using L1489PPC agreed within 25%. 

3. RESULTS 

3.1. L1495-B218 Filaments Seen in Am¬ 
monia 

3.1.1. Integrated intensity of NH 3 (1,1) & (%>%) 
and Definition of Subregions 

Figure [l] presents the integrated intensity of 
NH 3 (1,1) emission. The NH 3 (1,1) map clearly 
reveals dense cores associated with the L1495- 
B218 filaments. On the other hand, the filamen¬ 
tary structures between dense cores are barely de¬ 
tected, while those regions are very clearly visi¬ 
ble in an extinction map (Schmalzl et al. 2010) 
and 13 CO (Goldsmith et al. 2008; Pineda et al. 
2010b). This suggests that the density of fila¬ 
mentary structures between dense cores may not 
be high enough to form NH 3 efficiently (nn 2 < 
10 3 cm” 3 ; Flower et al. 2006) or that ammonia 
molecules exist but excitation conditions are insuf¬ 
ficient to observe the ( 1 , 1 ) emission at our average 
rms level of 120 rnK (Figure 2). The NH 3 (2,2) 
emission is detected only from dense cores and is 
much weaker than ( 1 , 1 ) emission, which implies 
that the kinetic temperature may be too low to 
excite to the J=2, K= 2 state except in dense gas 
of n Ha > 10 4 cm ” 3 (Ho & Townes 1983; Danby et 
al. 1988; Maret et al. 2009). 

We compare our NH 3 (1,1) emissions with 500 
/jm dust continuum emission seen by the SPIRE 
instrument of the Herschel Space Observatory in 
Figure [l] (Palmerim et al. 2013; Marsh et al. 


2014). The 500 p m dust continuum emission is 
thermal emission from dust grains and traces the 
column density structure of the filaments. The 
bright regions of 500 pm and NH 3 (1,1) emission 
show a good agreement with each other. The 500 
pm emission is more spatially extended than the 
NH 3 (1,1) emission. The lowest 500 pm intensity 
levels at which NH 3 (1,1) and (2,2) are detected 
at is 25 MJy/sr (Nh 2 = 5.6 x 10 21 cm ” 2 at 10 
K) and 177 MJy/sr (Nh 2 = 3.9 x 10 22 cm ” 2 at 
10K), respectively. This shows that NH 3 is a good 
tracer with which to identify intermediate to dense 
regions within filaments. 

In Figure [3j we zoom into each region and 
show the location of protostars identified using 
IRAS and the Spitzer Space Telescope (Rebull et 
al. 2010, 2011; Kryukova et al. 2012). We di¬ 
vide the map into subregions based on Barnard’s 
identification (Barnard 1927). We discuss each re¬ 
gion and determine if it is more evolved (having a 
Class I/flat spectra protostar associated with NH 3 
emission) or less evolved below. 

B7/L1495 contains multiple Class II and III 
protostars, two flat spectra protostars, and one 
Class I protostar, none of which is embedded 
within NH 3 core but there are nearby NH 3 cores. 
The three dense cores seen in NH 3 are mostly com¬ 
monly referred to as L1495A-S, L1495A-N, and 
L1495B (Lee et al. 2001). In B213 (also referred 
to as L1521D, Codella et al. 1997), there are two 
Class I protostars embedded in NH 3 cores in the 
eastern part of the region and another three Class 
I protostars, in the western part of the region. 
Therefore, B7 and B213 are the most evolved re¬ 
gions in the L1495-B218 filaments and B7 is likely 
to be more evolved than B213 (see Hacar & Tafalla 
et al. 2013). 

B218 has only one Flat spectra protostar em¬ 
bedded within a NH 3 core. Thus, we classify B218 
as an evolved region with star formation, although 
the NH 3 emission is dominated by two bright star¬ 
less cores. This region is also sometimes confus¬ 
ingly referred to as B217 (Benson & Myers 1989). 
We use the notation B218 in this paper since it is 
closer to the original Barnard catalog position of 
B218 (Kenyon et al. 2008). 

B10 also has only one protostar, which is Class 
II. In this work, we assume that B10 is a less 
evolved region because there is no other protostar 
except the one Class II protostar which is not di- 
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rectly associated with any dense cores in the re¬ 
gion. 

In contrast to the other regions, B211 and B216 
show weak NH 3 (1,1) emission and do not possess 
any known protostars. Thus, both regions are less 
evolved. B216 region is also named L1521B (Lee 
et al. 2001 ) and contains starless cores with large 
abundances of carbon chain molecules (Hirota et 
al. 2004). 

Hacar & Tafalla (2013) have also determined 
the relative ages of subregions using N 2 H + 1—0 
emission and the locations of protostars. Their 
estimates of relative age of subregions agrees well 
with ours except for B216. In their observation of 
N 2 H+ 1—0, they did not detect emission in B216. 
B216 has weaker NH 3 (1,1) emission and stronger 
CCS emission than B211 (Seo et al. in prep.) sug¬ 
gesting that B216 is a younger region. Thus, B211 
and B216 are both relatively young and inactive 
parts of the L1495-B218 filaments with B216 likely 
being the most newly condensed part of the fila¬ 
ments. 

3.1.2. Spectral line fitting 

To understand the physical properties of NH 3 
structures in the L1495-B218 filaments, we fit a 
simple NH 3 spectra model to the observed spectra 
and deduce physical quantities including kinetic 
temperature T^, excitation temperature T ex , total 
optical depth of NH 3 (1,1) emission n, dispersion 
of emission line a v , and LSR velocity vi sr (See Ap¬ 
pendix A for details of the fitting algorithm). The 
fitting algorithm includes only one velocity com¬ 
ponent, and this is satisfactory for all regions in 
the L1495-B218 filaments except B211 (see section 
3.2.2 for the LSR velocity of NH 3 ). We only fit¬ 
ted spectra where the detection of NH 3 (1,1) and 
( 2 , 2 ) emission is above 5<7T mb and 3ciT mb respec¬ 
tively, where crT mb is the baseline rms. 

3.1.3. Gas kinetic temperature and excitation 
temperature 

Maps of the gas kinetic temperature deduced 
from NH 3 (1,1) and (2,2) emission are presented 
in Figure [4j The gas kinetic temperature ranges 
from 8 K to 15 K in most of the region mapped. At 
the bright NH 3 peaks indicated by red arrows, the 
kinetic temperature decreases toward the core cen¬ 
ter. The drop in the kinetic temperature is about 


-1 to -3 K, or about 10% — 30% of the average 
kinetic temperature of surroundings. This tem¬ 
perature drop may ultimately break the dynami¬ 
cal stability of these dense cores and lead them to 
collapse (see Khesali et al. 2013). 

Figure [5] presents histograms of physical quan¬ 
tities to investigate statistical differences between 
subregions. In the first panel of Figure [5] we 
show histograms of the gas kinetic temperature. 
Medians of the gas kinetic temperature in sub- 
regions are all around 9.5 K with less than 1 K 
deviation (60% of all spectra have kinetic tem¬ 
perature within ±1 K of 9.5 K.). The two less 
evolved regions, B10 and B216, have lower kinetic 
temperature than the other more evolved regions 
(we exclude B211 from discussion because it has 
multiple velocity components and our NH 3 fit¬ 
ting model works only for a single velocity compo¬ 
nent.). The histograms of the gas kinetic temper¬ 
ature demonstrate a tendency that a less evolved 
region is cooler than a more evolved region and 
that a denser region is cooler than diffuse region 
if there are no embedded protostars; however, the 
gas kinetic temperature differences are subtle. 

Histograms of NH 3 excitation temperature are 
shown in the second panel of Figure [5] Median 
values of the excitation temperature are mostly 
around 5 K, which is only half of the median val¬ 
ues of the kinetic temperature. This demonstrates 
that NH 3 in all regions is subthermally excited. 

3.1.4- LSR velocity 

NH 3 lines at 39 NH 3 integrated intensity peaks 
in the L1495-B218 filaments are shown in Figures 
[6]and[7](for peak locations, see §3.2.1). The range 
of NH 3 vi sr is from 5.0 km s -1 to 7.2 km s” 1 , 
which is the same range seen in the C 18 0 J = 
1—0 spectra of Hacar et al. (2013). NH 3 struc¬ 
tures are confined to a single velocity component 
except in B211. Comparing NH 3 to C 18 0 (Hacar 
et al. 2013) and 13 CO (Goldsmith et al. 2008), 
we found that all NH 3 peaks have corresponding 
velocity components in 13 CO, while only 28 out of 
39 NH 3 peaks directly associated with the same 
velocity components in C 18 0. Some NH 3 peaks 
where NH 3 integrated intensity is relatively high 
(>40 K km s _1 , i.e. four core-like structures neigh¬ 
boring in a row in B213) do not have correspond¬ 
ing velocity components in C 18 0 and 13 CO emis¬ 
sion at those places is also weak. We also found a 
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few NH 3 peaks with both NH 3 and C 18 0 emission 
but discrepancies in their LSR velocities, which 
may be a sign of chemical differentiation. 

The LSR velocity structure of the L1495-B218 
filaments was extensively studied by Hacar et 
al. (2013) using C 18 0 J= 1—0. They identi¬ 
fied 35 filaments using a “friend of friends” al¬ 
gorithm (FIVe; Hacar et al. 2013) and for each 
filament they deduced physical properties includ¬ 
ing mass, dispersion of i>i sri nonthermal motion, 
mean kinetic temperature. (See Table 3 in Hacar 
et al. 2013). Using dense cores identified from 
their N 2 H + J = 1—0 observations, they suggested 
only 7 out of 35 filaments are “fertile” meaning 
they contain a N 2 H + -briglit dense core. How¬ 
ever, N 2 H+ traces gas with density higher than 
10 4 cm -3 , which is the central density of rela¬ 
tively evolved dense cores (Lada et al. 2008). It 
is likely that dense cores at early stages are not 
fully revealed by N 2 H + . Ammonia molecules have 
an intermediate effective excitation density that 
is between C 18 0 J = 1—0 and N 2 H + J = 1—0; 
therefore, our ammonia map may probe the ve¬ 
locity structure of both developing dense cores as 
well as later stages. 

Matching the LSR velocity of NH 3 (1,1) and 
(2,2) structures to 35 filaments of Hacar et al. 
(2013), we found NH 3 (1,1) emission in filaments 
6 , 8 , 10, 11, 15, 20, 28, 32 (Table 3 of Hacar et 
al. 2013). Hacar et al. (2013) suggested that 
filaments 6 , 10, 11, 15, 20, 32 and 34 are fertile 
( i.e. having dense cores in their N 2 H + observa¬ 
tion). Filament 34 is not within the field of view 
of our observation, so it is excluded in this study. 
Our NH 3 observation reveals that two more fila¬ 
ments, 8 and 28, also have NH 3 emission. Hacar 
et al. (2013) argued that a filament tends to be 
sterile when the line density of filament is less than 
15 M 0 pc -1 , which is the critical density for an 
isothermal cylinder at 10 K (Stodolkiewicz 1963; 
Ostriker 1964). Filament 28 has a line density of 
24.1 M 0 pc -1 (Hacar et al. 2013) but there was no 
detection of dense cores in N 2 H + J = 1—0. With 
our NH 3 observations, we found that there are at 
least two dense cores in filament 28. On the other 
hand, filament 8 has a line density of 9.3 M 0 pc -1 
(Hacar et al. 2013), which is lower than the criti¬ 
cal line density, but we see that there is a relatively 
bright NH 3 core. Since Hacar et al. (2013) esti¬ 
mated the masses of filaments from C 18 0 J = 1—0, 


it may be possible that they underestimated the 
mass due to CO depletion. The filament mass es¬ 
timated from dust continuum may be larger than 
mass estimated from C ls O J = 1—0, but the fil¬ 
aments cannot be distinguished along the line of 
sight, which will result in an overestimation of fil¬ 
ament masses. From the fact that filament 8 is 
relatively bright in NH 3 (1,1), it is likely that the 
line density of 9.3 M 0 pc -1 is the minimum line 
density of filament 8 and the true line density may 
be higher than the critical value. 


3.1.5. Velocity dispersion 

Velocity dispersions range from 0.05 to 0.25 km 
s _1 except in B211 (B211 has multiple velocity 
components which are not pursued in this study). 
Overall, we found that velocity dispersions, which 
include both thermal and nonthermal components, 
tend to be wider in the less evolved regions than 
those in the more evolved regions. Histograms of 
the velocity dispersions (Figure [ 5 ]) show this trend 
clearly. The median values of velocity dispersions 
for the least evolved regions, B216, is 0.118“ 4 , 
while the most evolved regions, B7, has a median 
dispersion equal to 0.104 km s _1 (the thermal ve¬ 
locity dispersion of NH 3 at 10 K is 0.069 km s _1 ). 

Using accurate gas kinetic temperatures de¬ 
duced from NH 3 (1,1) and (2,2), we estimate the 
velocity dispersions of nonthermal kinetic motions 
er nt (Figure [8]) as follows: 


&nt — 



kT k 

Mnh 3 mh ’ 


(1) 


where k is Boltzmann constant, /tnh 3 is the molec¬ 
ular weight of ammonia (17 a.m.u.), and mh is the 
mass of a hydrogen atom. We find most regions 
of the filaments are supported mainly by thermal 
pressure (the fourth panel in Figure [5]) while some 
edges of filaments are dominated by nonthermal 
support (cr„t > ^thermal = 0.188 km s _1 at T k = 
10 K and p = 2.33). The median ratio of ther¬ 
mal support to nonthermal support ( R p = c 2 /a 2 t ) 
is around five. In addition, the median ratios of 
thermal support to nonthermal support tends to 
be larger for more evolved regions (median R p = 
6.3 in more evolved region B7, while R p = 3.7 
in less evolved region B216). This suggests that 
nonthermal motions are considerably dissipated in 
more evolved regions and making it easier to ini- 
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tiate star formation, which has been seen toward 
dense cores in other molecular clouds (e.g. Good¬ 
man et al. 1998; Pineda et al. 2010a; Tatematsu 
et al. 2014). 

3.2. Dense NH 3 Structures Within the 
L1495-B218 Filaments 

3.2.1. Identifying dense structures in the L1495- 
B218 filaments 

The L1495-B218 filaments have many dense 
cores in different evolutionary stages of star for¬ 
mation, ranging from starless to protostellar cores. 
We use our ammonia observations to identify 
dense cores in the L1495-B218 filaments using the 
Cardiff Source-finding AlgoRithm (CSAR, Kirk et 
al. 2013), which uses both a seeded watershed and 
dendrogram algorithm to find clumpy structures 
and their spatial relationships. We used the inte¬ 
grated intensity of the central hyperfine group (the 
strongest hyperfine group of nitrogen splitting) of 
NH 3 (1,1) to identify dense structures because the 
central hyperfine group gives the highest signal- 
to-noise ratio compared to other four hyperfine 
groups. Parameters in the CSAR algorithm are 
set to find structures above a 7.5-cr integrated in¬ 
tensity level in the (1,1) transition (~3cr in (2,2) 
transition) and to identify a clump with a size 
larger than 9 m b/2 ~ 15" (fullbeam option) and a 

2 - cr enhancement in intensity from its background. 
We also used the removeedge option in order not 
to have overextended edges for a source. In this 
study, we use 2 -cr as an intensity interval to clip 
intensity peaks instead of using typical value of 

3- cr because with 3-cr intensity interval CSAR al¬ 
gorithm tends to combine two peaks as one elon¬ 
gated clumpy structure. These peaks typically lie 
on top of nested sources which are many a brighter 
than their background. The CSAR algorithm also 
had trouble in excluding intensity spikes at edges 
of maps which have a high noise level ( ar k > 0.3 
K). In order to exclude the edges of NH 3 maps, we 
made 12 patches and applied the CSAR algorithm 
to each patch. 

The main reason we use NH 3 (1,1) instead of 
dust continuum is because NH 3 emission acts as a 
volume density filter. From the excitation curve of 
NH 3 (1,1), the effective excitation density varies 
slowly over the range of gas kinetic temperature 
from 7 K to 15 K. The excitation temperature of 


dense structures identified by the CSAR algorithm 
is typically above 3.2 K, and the effective excita¬ 
tion density for the range of observed gas kinetic 
temperature and NH 3 abundance varies only from 
7xl0 2 cm -3 to 1 xlO 3 cm -3 . This suggests that 
NH 3 (1,1) emission naturally filters out regions 
with a density lower than ~ 7 x 10 2 cm -3 . This 
may make NH 3 structures appear smaller than 
corresponding dust continuum structures, but the 
NH 3 sources identified by the CSAR algorithm are 
true structures in volume density whereas analy¬ 
sis in dust continuum sees only an enhancement 
in total column density. 

From the integrated intensity of the central hy¬ 
perfine group of NH 3 (1,1), we found 39 leaves 
(a leaf: a structure containing an intensity peak) 
and 16 branches (a branch: a nested group of 
leaves). Figure [9] sh ows the locations of NH 3 
leaves, and FigurejlOl shows hierarchy among NH 3 
leaves and branches in a dendrogram. The den¬ 
drogram has no lowest common branch because 
NH 3 (1,1) structures are highly fragmented due 
to the volume density and chemical filtering. In 
this study we refer to “active” leaves as those 
which have a Class I or Flat spectrum protostar 
within their lowest-level branches in the dendro¬ 
gram (a.k.a. root) or within a distance of three 
times their radius. Note that this definition means 
that active leaves could be either starless or pro¬ 
tostellar. Only 4 leaves have embedded protostars 
and 17 leaves are classified as active. Hacar et 
al. (2013) found 19 dense cores in N 2 H+ 1—0, 
while CSAR finds 39 leaves in NH 3 (1,1). This 
difference may be because NH 3 is an intermediate 
density gas tracer that reveals more cores than the 
denser gas tracer N 2 H + and our NH 3 observation 
($mb ~ 31") have better angular resolution than 
the N 2 H + observation (0 m b ~ 50"). 

For the comparison between dust and NH 3 , we 
also identified dust cores using CSAR. From the 
500 /inr dust continuum data of the Herschel Space 
Observatory (0 m b ~ 36"), we found 51 dust leaves. 
Figure [ll] shows the locations and boundaries of 
dust leaves. 32 dust leaves coincide with 38 NH 3 
leaves, while 19 dust leaves do not agree with any 
NH 3 leaf. Nine of the dust leaves contain at least 
two NH 3 leaves. This suggests that there may be 
multiple dense cores or fragments within a core 
since NH 3 is sensitive to volume density whereas 
the dust continuum is optically thin and depends 
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on the total column density. 

In order to compare how well dust continuum 
and NH 3 agree with each other, we investigate the 
separation of the dust peak position and the NH 3 
peak position in a dust leaf associated with a NH 3 
leaf. Most of NH 3 peak positions show agreement 
with corresponding dust peak positions with a me¬ 
dian separation of 18".7, which is slightly larger 
than half the FWHM beam size. However, there 
are four dust leaves having a separation between 
dust and NH 3 peaks larger than 30" (the maxi¬ 
mum separation is 66".2), and three out of four 
are associated with a nearby protostar (Class I or 
Flat spectrum). This suggests a large separation 
may be due to chemical effect produced by the ra¬ 
diation from a protostar. The fact that most of the 
NH 3 sources agree with the dust sources demon¬ 
strates that NH 3 sources identified by CSAR may 
represent dense structures within the L1495-B218 
filaments. In the following subsections all of the 
analysis are done with NH 3 sources determined 
from our NH 3 data. 


3.2.2. Size and mass 


The physical properties of NH 3 leaves and 
branches are summarized in Tables 0 0 and 0 
The quoted kinetic temperature and velocity dis¬ 
persion are the median values estimated from NH 3 
fittings within NH 3 leaves and branches except for 
NH 3 leaf 11. This has a peak NH 3 (1,1) line inten¬ 
sity weaker than 5cr except at its peak location, so 
we fit the NH 3 line at the peak intensity location. 

The column density of molecular hydrogen 
within each structure is estimated from the 500 
/jm continuum of Herschel Space Observatory us¬ 
ing the following equation. 


^h 2 




^gas 


(^dust)^i/^ap ^Idust 


( 2 ) 


where 7 Vh 2 is column density of H 2 molecules, S v is 
observed flux at frequency zq p, is the mean molec¬ 
ular weight, mh is the mass of hydrogen atom, B v 
is black body radiation at frequency zq T<j ust is 
temperature of a dust grain, is opacity of dust 
grain at frequency zq P is the solid angle used in 
the observation, and mdust / m g!is is dust-to-gas ra¬ 
tio, which we set to 0.01. In this work, we assume 
Tdust = Tj t, where T% is the kinetic temperature 
deduced from our ammonia observations, p is 2.8. 


We use k u of 5.04 cm 2 g -1 (OH5), which is the 
column 5 of “MNR with thin ice mantle” model 
from Ossenkopf & Henning (1994). The peak H 2 
column density of NH 3 sources in the L1495-B218 
filaments ranges from 6.0 x 10 21 to 4.3 x 10 22 
cm~ 2 , which is similar to the typical dense core 
column density in other molecular clouds (e.g. Ji- 
jina et al. 1999; Rosolowsky et al. 2008). 


We define the size of a NH 3 source as the radius 
of a circle that has an area equal to that of a NH 3 
CSAR leaf contour. In the L1495-B218 filaments, 
the FWHM from Gaussian fitting may not be a 
useful representation of the size of a dense core 
because the FWHM contour of a dense core often 
overlaps with the FWHM contour of its neighbor¬ 
ing dense core. Since CSAR separates neighboring 
sources through hierarchical structure, we use the 
area identified by CSAR to estimate the size of a 
structure. Structures were found on scales ranging 
from 2x 10” 2 pc (4375 AU) to 0.1 pc (20000 AU). 

The masses of NH 3 sources are estimated by 
summing up the column density deduced from 500 
pan continuum within each NH 3 leaf and branch 
(Table 0 and black circles in Figure 10). This is 
a typical method for estimating mass of a dense 
source since three dimensional morphology is usu¬ 
ally unknown. If three dimensional morphology is 
known or assumed, mass of a local background 
may be subtracted, which will result a smaller 
dense source mass than the dense source mass esti¬ 
mated by simply summing up the column density. 
The resulting masses of NH 3 leaves range from 
0.05 Mq to 4.3 M q , which are similar to typical 
dense core masses observed in other studies (0.1-10 
M 0 ; Motte et al. 1998, 2001; Schmalzl et al. 2010; 
Onishi et al. 1996, 2002; Kirk et al. 2013). The 
masses of NH 3 branches range from 0.78 Mq to 9.5 
Mq. In addition, we also compared three differ¬ 
ent methods of determining masses of NH 3 leaves 
in Table 0 and Figure 12 (mass estimated from 
dust continuum, virial mass shown in §4.1, and 
mass estimated from NH 3 column density with 
the fractional abundance of NH 3 relative to hydro¬ 
gen molecules being 1.5 x 10 -9 ). The masses are 
strongly correlated with each other with a Spear- 
mann’s rank correlation coefficient of p = 0.93. 
This strong linear correlation suggests that the as¬ 
sumptions (i.e. opacity and temperature) used to 
derive mass from the dust emission do not have 
strong systematic variations from the lowest mass 





to highest mass leaves. It also suggests that the 
properties derived from NH 3 sources identified by 
CSAR are representative of the properties of dense 
regions in the L1495-B218 filaments, regardless of 
whether those regions are traced by dust emission 
or NH 3 . 

We also compared nonthermal support to ther¬ 
mal support within NH 3 leaves identified by CSAR 
in Figure[l2j All of the NH 3 sources are dominated 
by thermal support, with an average value of R p 
= Cs/cr,j t equal to 4.9, and there is no dependency 
on R p vs. mass. This shows that nonthermal mo¬ 
tions are efficiently dissipated in dense regions and 
making it easier to initiate star formation. 

3.2.3. Shape and orientation 

The two dimensional morphology of the NH 3 
sources is estimated using the moment of inte¬ 
grated intensity of NH 3 (1,1) with the center at 
the location of the intensity peak within a leaf or 
a branch. From the moment of integrated inten¬ 
sity we derive the principal axes, the axis ratios, 
and the angles with respect to galactic longitude, 
which are also summarized in Table [3 The dis¬ 
tribution of angles of NH 3 leaves with respect to 
their lowest-level branches in dendrogram is shown 
in the third panel of Figure [13| The minimum and 
the maximum phase angles of NH 3 leaves with re¬ 
spect to their lowest-level branches are -49° and 
53°, respectively. Most of NH 3 leaves are well 
aligned with their lowest-level branches within a 
range of ±40°. Since a lowest-level branch cor¬ 
responds to a filament, this suggests that NH 3 
leaves are well aligned with filaments and were 
likely formed within the filaments. 

The axis ratio of the NH 3 leaves is estimated 
as short axis/long axis (Table [3]). The distribu¬ 
tion of axis ratios is shown in the fourth panel of 
Figure [i~3} The axis ratios range from 0.2 to 1.0 
and are mostly around 0.5, which is close to the 
average value obtained from dense core observa¬ 
tions by Myers et al. (1991). In order to deduce 
the three dimensional morphology of NH 3 leaves, 
we compared the observed axis ratio distribution 
to the probability distribution of apparent axis ra¬ 
tios of spheroids using a Monte Carlo method (see 
Appendix B for more details). We found the best- 
fit axis ratio distribution when the spheroids are 
prolate and the orientation angles (the angle be¬ 
tween the equatorial planes of spheroids and the 


observer’s line of sight) range from -80° to 80° 
(thick solid line in the fourth panel of Figure [l3|- 
The parameters of the prolate spheroids of the 
best-fit model are a mean axis ratio (short axis/ 
long axis) of 0.40 and a dispersion of 0.11 with \ 2 
= 5.92 . We could not fit the distribution model 
using oblate spheroids within 3 -a uncertainty (x 2 
= 9) unless we limit the orientation angles much 
narrower than the observed distribution of the an¬ 
gles between the axes of NH 3 leaves and the axes 
of filaments (third panel of Figure 131. Thus, the 
NH 3 leaves are more likely prolate spheroids than 
oblate spheroids. 


4. Discussion 

The physical and dynamical properties of dense 
structures depends on the history of their forma¬ 
tion. Theoretical studies of molecular filament 
and dense core formation in a turbulent envi¬ 
ronment (the gravoturbulence fragmentation sce¬ 
nario) shows that a molecular cloud fragments and 
develops local density enhancements by colliding 
flows of turbulence and the density enhancements 
are grown by self-gravity to evolve into filaments 
and dense cores. Dense structures formed follow¬ 
ing this scenario tends to be pressure-confined in 
their early evolution stage and proceed to a gravi¬ 
tationally bound state. On the other hand, dense 
cores formed in a quiescent, magnetized environ¬ 
ment are typically gravitationally bound through 
their evolution since fragmentation typically takes 
place when a structure loses support against grav¬ 
ity. In addition, three dimensional morphologies 
and alignments of dense structures with magnetic 
fields also depend on the history of their forma¬ 
tion. In a quiescent environment, magnetic fields 
plays an important role in shaping and aligning 
dense structures, while the magnetic fields may 
not play a major role in determining the prop¬ 
erties of dense structures in a turbulent environ¬ 
ment. In this section, we discuss the dynamical 
status, three dimensional morphologies, and sta¬ 
tistical properties of dense structures in the L1495- 
B218 filaments as well as uncertainties in our anal¬ 
ysis. 

4.1. Virialization 

We investigate the dynamical states of NH 3 
leaves using the virial theorem. An equation of 


9 



the virial theorem without magnetic contributions 
is 


~ = W + 2(T + H)~ Jpx-MS (3) 

where I is the moment of inertia, W is the gravi¬ 
tational energy, T is the kinetic energy of systemic 
motions, II is the internal energy of thermal and 
nonthermal motions, and the integral term is the 
work done by the external pressure p. S is the core 
surface and fi is the unit vector perpendicular to 
the core surface. We omitted the terms related to 
magnetic fields because there is no measurement 
of magnetic fields within NH 3 leaves. There are 
studies using the virial theorem including all terms 
except magnetic fields (e.g. Hunter 1979; Seo et al. 
2013) but only the gravitational and the internal 
energies are conventionally used for studying dy¬ 
namical stability of a dense core in observational 
studies because it is hard to estimate other terms. 
We also start with only the gravitational and the 
internal energies in estimating the virial mass and 
further check whether or not the pressure term is 
important in confining dense cores. We assume a 
time-independent steady state so that cPl/dt 2 = 
0 . 

We assumed a NH 3 leaf as a homoeoidal el¬ 
lipsoid with density falling as p~r~ 3 (Bertoldi & 
McKee 1992). The gravitational energy of a ho¬ 
moeoidal ellipsoid is 


TTr 3 1 — f GM 2 arcsine 

W = -A- 

5i-M R e 

5 


(4) 


where e is the eccentricity of homoeoidal ellipsoid, 
and R is the semi-major axis. The internal energy 
is estimated by 



(5) 


where d v is a representative value of velocity dis¬ 
persion including thermal and nonthermal mo¬ 
tions within a NH 3 leaf. It can be either the av¬ 
erage, median, or mass-weighted values. In this 
study, we use median values of a v within NH 3 
leaves. Finally, the virial mass when W = 2 U 
is given as 


Mvir = 5 


1 - f Rdl 


1 — i. G arcsin e 
3 


( 6 ) 


We first estimated a virial mass with p ~ r ~ 2 
(/3 = 2), which is listed in Table [ 2 ] The ec¬ 
centricities of leaves are estimated from the axis 
ratios of NH 3 leaves which are given in Table [3] 
The virial masses of NH 3 leaves are shown as red 
circles in Figure 10 The virial mass is calcu¬ 


lated for only the NH 3 leaves because the NH 3 
branches usually have either very elongated or ir¬ 
regular shapes, which makes it difficult to cal¬ 
culate the virial mass with conventional assump¬ 
tions. We find that 30 leaves out of 39 leaves have 
observed masses smaller than their corresponding 
virial masses (the first panel in Figure |~i~3| ), which 
suggests that most of leaves are gravitationally un¬ 
bound (the median ratio is about two). 7 out of 9 
gravitationally bound leaves are active (red), while 
only two out of 22 inactive (blue) leaves are grav¬ 
itationally bound. A virial mass of a homoeoidal 
ellipsoid with is 50% larger than the cor¬ 

responding virial mass of a homoeoidal ellipsoid 


with For p~r _i , there are only two grav¬ 

itationally bound leaves, both of which are active 
leaves. 


In order to assess the validity of the virial anal¬ 
ysis, we must first carefully analyze the uncertain¬ 
ties. One uncertainty may come from a choice of 
d v . In this study, we chose median values of a v 
within NH 3 leaves but one may use a mean of a v or 
a mass-weighted mean of cr v (mass weighted mean 
is d 2 = / NYL 2 pm h o 2 dAI J N H2 pm h dA ) within a 
NH 3 leaf. We found that the mean values of a v 
are at most 10% larger (typically 5%) than the 
corresponding median values, which may result in 
at most 20 % increase of corresponding virial mass. 
The mean values may be larger than the median 
values because a v tends to be lower at the center 
of a NH 3 leaf and considerably increase toward 
edges of the NH 3 leaf. Mass weighted means on 
the other hand are in the range from 95% to 105% 
of the corresponding median a v values. Thus, un¬ 
certainties originating from choices of d v do not 
strongly affect the conclusion of virial mass analy¬ 
sis since ratios of viral masses to observed masses 
are quite large with respect to these uncertainties. 

The largest uncertainties in estimating the ob¬ 
served mass from 500 /jm continuum data are 
mainly due to assuming Td = Tk and uncertainty 
in dust opacity n. Comparing the dust tempera¬ 
ture, which is estimated by SED fitting of 70, 160, 
250, 350, and 500 /mi dust continuums (Palmeirim 
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et al. 2013), to the gas kinetic temperature de¬ 
duced from NH 3 (1,1) and (2,2), we found that 
the dust temperature is typically 3 — 5 K higher 
than the gas kinetic temperature. If the dust tem¬ 
perature is 5 K higher than the gas kinetic tem¬ 
perature, the mass estimated from the dust tem¬ 
perature is 35% less than the mass estimated with 
the Td = Tfe assumption and only four active NH 3 
leaves are gravitationally bound. 

Another uncertainty is from the dust-to-gas ra¬ 
tio. In this study, we use 0.01, which is a standard 
value assumed in a molecular cloud (e.g. Beckwith 
et al. 1990; Lombardi et al. 2006; Schmalzl et al. 
2010; Kirk et al. 2013). However, detailed ob¬ 
servations of interstellar dust grains in the diffuse 
ISM indicates that the dust-to-gas ratio is close to 
0.0064 (Draine et al. 2011 and reference therein). 
This increases observed mass of the leaves about 
55% from the current values. In this case, 20 out of 
39 leaves are gravitationally bound and 19 leaves 
are still gravitationally unbound. 

Not knowing three dimensional morphologies of 
leaves also brings uncertainty in estimating masses 
of leaves. We estimate the mass of a leaf by sim¬ 
ply summing up the column mass of pixels that 
belongs to the leaf, which corresponds to a cylin¬ 
drical morphology rather than a spheroidal mor¬ 
phology. This has been a conventional method 
to estimate the masses of dense structures since 
we do not know three dimensional morphologies 
(e.g. Kirk et al. 2013). In order to check the un¬ 
certainty from assuming different morphologies in 
estimating the mass of a dense structure, we es¬ 
timated the masses of leaves by summing up the 
column density of pixels within leaves after sub¬ 
tracting the average column density of local back¬ 
grounds/branches. We found that the masses of 
leaves are typically 20 — 50% less than the origi¬ 
nal values if we subtract a local background from 
the leaves in estimating leaf mass and that only 
two our of 39 leaves are gravitationally bound. 

The typical uncertainty in the dust opacity is 
a factor of few (e.g. Shirley et al. 2011) and this 
results in quite a large uncertainty in estimating 
mass since M oc 1/ft. We use OH5 dust opacity but 
this is not a direct measurement of dust opacity 
from molecular clouds but a simulated dust opac¬ 
ity based on a coagulated dust model. Shirley et 
al. ( 2011 ) made a direct measurement of the dust 
opacity at 450 /im and 850 pm in the outer enve¬ 


lope of a low-mass Class 0 core, B335, through a 
comparison between infrared extinction and sub¬ 
millimeter dust continuum. We compared the dust 
opacity of Shirley et al. (2011) to OH5 dust opac¬ 
ity at 450 /jm and applied the same uncertainty 
range of dust opacity to OH5 dust opacity at 500 
/jm (ft = S.Od'L'f^g cm 2 g -1 ). The uncertainty in 
the observed mass ranges from 57% lower to 28% 
higher with respect to the observed mass given in 
Table [2] If we take the lower limit of uncertainty 
in the observed mass, two active NH 3 leaves are 
gravitationally bound, while if we take the upper 
limit of uncertainty in the observed mass 15 out 
of 39 are gravitationally bound. 

Thus, uncertainties in a v , Td, md/m gi geome¬ 
try, backgrounds, and dust opacity are unlikely to 
change our conclusion that most leaves are grav¬ 
itationally unbound. This result suggests that a 
dense core may firstly form as a gravitationally un¬ 
bound structure, evolve to a gravitationally bound 
core, and then undergo collapse to form a proto¬ 
star. Indeed, 10 out of 15 bound leaves in the 
extreme opacity limit are in more evolved regions. 

4.2. Pressure-confinement 

External pressure may be a source to confine 
gravitationally unbound structure. Dib et al. 
(2007) carried out a detailed virial analysis for 
dense cores forming in a turbulent environment 
using a numerical method. They showed that the 
work done by the external pressure in the virial 
theorem is comparable to the other terms such as 
internal energy and gravitational energy (Dib et 
al. 2007). Since most of the NH 3 leaves are not 
gravitationally bound in our simple virial analysis, 
we checked whether or not work done by the ex¬ 
ternal pressure is comparable to the internal and 
gravitational energy. The work done by external 
pressure is estimated as follows 

px ■ MS = 4 — mm{N H2 idm h al} c A, (7) 

where A is the projected area of NH 3 leaf. 
Nn 2 pm.hO'v is the integration of pressure along 
the line of sight, where Nh 2 is column density 
estimated from 500 pm and a v is velocity disper¬ 
sion of molecular gas with p = 2.33. We chose 
the minimum value at the circumference of NH 3 
leaf, min{A^H 2 ^ni^cr 2 }c, because it corresponds 
to the minimum external pressure exerted on a 
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NH 3 leaf. The coefficient of 4/(3 a/ 1 — e 2 ) is for 
spheroid geometry when one of the principal axes 
is aligned with the observer’s line of sight. We 
estimate the work done by external pressure using 
min{iVH 2 p.mfc(7^}c because the external pressure 
is not well known. The estimated work done by 
external pressure ranges from 5 x 10 40 erg to 1 
x 10 42 erg. In a majority of leaves, gravitational 
energies are larger than the corresponding work 
done by external pressure except for 12 leaves 
(the second panel of Figure 13), which suggests 
that those 12 leaves are pressure-confined. Active 
leaves tend to have a lower average ratio than that 
of inactive leaves, which indicates the gravity be¬ 
comes more important than external pressure as a 
dense core evolves. Moreover, in 9 out of 12 leaves, 
the ratios of the work done by external pressure 
to the gravitational energy are larger than that 
of the critical Bonnor-Ebert sphere (black dashed 
vertical line). Thus, these results suggest that the 
youngest condensations are pressure confined. 


A source of the external pressure may be inves¬ 
tigated by comparing the internal pressure of the 
NH 3 leaves and the pressure of the surrounding 
media. The mean internal pressure of NH 3 leaves 
is < p/k > = 2 x 10 6 K cm -3 , where p is roughly 
estimated as 


P = 


Mol 

V 


( 8 ) 


where V is the leaf volume and a v is the median 
value of a v within a NH 3 leaf. We assumed that 
leaf volume is the same as the volume of a pro¬ 
late spheroid having the observed axis ratio (see 
§4.2.3 for shapes of NH 3 leaves). If we assume the 
density structure of a leaf resembles the critical 
Bonnor-Ebert sphere, the surface pressure is about 
40% of the mean internal pressure, < p su rface/fc > 
= 8 x 10 5 K cm -3 . This is significantly higher 
than the gas pressure of the interstellar medium 
(< PiSM/k > « 10 4 K cm -3 , Bertoldi & McKee 
1992; Draine 2011) and still a bit higher than the 
interstellar turbulence ram pressure, < p lam /k > 
a; 5 x 10 4 K cm -3 (Lada et al. 2008). 

Another source of external pressure may be the 
weight of the filament material outside of NH 3 
leaves since the self-gravity of the filaments pro¬ 
duces a gravitational potential well and the weight 
of filaments material outside of NH 3 leaves ex¬ 
erts a confining pressure on the NH 3 leaves. We 


estimated the pressure due to the weight of the 
filament material outside of the NH 3 leaves as 
(Bertoldi & McKee 1992; Lada et al. 2008): 

Pfilaments = " (9) 

where E is the mean mass surface density of mate¬ 
rial outside of leaves, and a is a correction factor 
for the non-spherical geometry of the filaments. 
The values of a range from 1.0 to 3.3 depending 
on the eccentricities and density structures of the 
filaments (for a = 3.3, we took the maximum axis 
ratio of the NH 3 branches, max(Z/R) ~ 10, and 
a density structure of p ~ r~ 2 ). The mean of 
the pressure of the filaments due to their weight 
is < Pfliament/fc > = 6 x 10 5 K cm -3 , which is 
a bit smaller than the surface pressure estimated 
from Equation ([ 8 ]) but within a factor of two. This 
suggests that the weight of the filaments is one 
of main sources that confines gravitationally un¬ 
bound leaves. 

Ram pressure due to the inflow of materi¬ 
als onto the filaments and dense cores from the 
molecular cloud may be another source of confin¬ 
ing dense structures. Theoretical studies demon¬ 
strated that ram pressure due to mass accretion 
to dense structures may provide considerable con¬ 
fining pressure (Heitsch 2013a, 2013b) if a con¬ 
verging flow or an inflow is relatively isotropic (a 
non-isotropic converging flow may disrupt a dense 
structure, Did et al. 2007). 12 CO and 13 CO ob¬ 
servations of the L1495-B218 filaments show that 
there is a difference of 2 km s _1 in the LSR veloci¬ 
ties from the north-east to the south-west of B211 
(Goldsmith et al. 2008; Palmeirim et al. 2013). 
This is similar to an infall velocity predicted in 
a similarity solution of gravitational infall onto 
a cylinder (Kawachi & Hanawa 1998). If we ac¬ 
cept that the difference in the LSR velocity is due 
to relatively isotropic inflow of materials from the 
molecular cloud to the filaments, the ram pressue 
due to an inflow motion is 

Paccretion = Pp(R)3V , (10) 

where p p is the density of the best-fit Plummer 
model to the filament (Palmeirim et al. 2013), 
and R is the outer radius of the filament which is 
0.4 pc for the B211-B213 regions. SV is the dif¬ 
ference in the LSR velocity across B211 and B213 
regions, which is about 2 km s _1 . The ram pres¬ 
sure due to the accretion of materials from the 
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molecular cloud to filaments are < pmament/fc > 
~ 4.5 x 10 5 K cm -3 , which is smaller than the 
surface pressure estimated from Equation (| 8 j) but 
still within a factor of two. Since there is projec¬ 
tion effect in measuring the true inflow velocity, 
this may be a lower limit of ram pressure due to 
an inflow motion. Thus, ram pressure due to in¬ 
flow of materials from the molecular cloud to the 
filaments may also be one of main sources of pres¬ 
sure that confine gravitationally unbound leaves. 

4.3. Aligned prolate dense structures 

The formation of a filament or a dense core 
through the gravoturbulence fragmentation sce¬ 
nario predicts that a dense structure may not be 
strongly aligned with magnetic fields (Nakamura 
& Li 2008) unless the magnetic field pressure is 
considerably stronger than the ram pressure of 
turbulence. On the other hand, formation through 
a quasi-static gravitational contraction may result 
in a strong alignment of the dense core with mag¬ 
netic fields because contraction occurs more effi¬ 
ciently along magnetic field lines. Although mag¬ 
netic field directions are not measured within the 
NH 3 sources, magnetic field directions on filament 
scales (a few 0.1 pc) have been measured through 
dust polarization in optical and infrared wave¬ 
lengths along the L1495-B218 filaments (Chapman 
et al. 2011 and references therein). In B7 and 
BIO the magnetic fields are relatively parallel to 
the long axes of a filament. From B213 to B216 
the magnetic fields are perpendicular to the long 
axes of filaments. Dust polarization in B218 is not 
measured. Since we do not know the magnetic 
field directions within dense cores and magnetic 
field directions are aligned in only some regions, it 
is difficult to ascertain the role of magnetic field 
in aligning dense cores. 

In the gravoturbulence fragmentation scenario, 
a dense core forms as an oblate spheroid and 
quickly proceeds to a prolate spheroid or triax- 
ial shape (e.g. Gong & Ostriker 2011; Kainu- 
lainen et al. 2014). The ratio of oblate dense 
cores to prolate dense cores depends on the Mach 
number of turbulence and the evolution time. A 
system with either a higher Mach number or a 
younger age shows a higher proportion of oblate 
dense cores. Since the NH 3 leaves are likely to 
be prolate (see §3.2.3), the L1495-B218 filaments 
may have been formed in a low Mach number tur¬ 


bulence (M < 2; Gong & Ostriker 2011) if the 
filaments are formed through the gravoturbulence 
fragmentation scenario. The velocity dispersions 
of filaments in the Taurus estimated from 13 CO 
observations are indeed low Mach number disper¬ 
sions that range from 0.15 km s _1 to 0.65 km s -1 
(Panopoulou et al. 2014), corresponding to Mach 
numbers of 0.65 — 3.2, if we take the kinetic tem¬ 
perature of 15 K (Goldsmith et al. 2008). Alterna¬ 
tively, magnetic fields may also make a dense core 
to have a non-spherical shape because a dense core 
tends to contract along the magnetic field lines 
whereas a direction perpendicular to the field is 
supported by magnetic pressure. If magnetic pres¬ 
sure is the main force making a dense core’s non- 
spherical shape, the dominant shape is an oblate 
spheroid (Lizano & Shu 1989). Since most of the 
NH 3 leaves are likely prolate spheroids, magnetic 
pressure does not seem to be the main cause of 
shaping a dense core into an elongated morphol¬ 
ogy, but gravoturbulence fragmentation may have 
played important roles in shaping the NH 3 leaves. 

4.4. Comparing Size, Mass, and Velocity 
Dispersion 

Larson proposed that there is an empirical rela¬ 
tion between line width and size of a structure in a 
cloud scale (Larson 1981). To see whether there is 
such a relation in a sub-filament scale, we compare 
the mean velocity dispersion of NH 3 sources with 
the NH 3 source size in Figure[l4j The velocity dis¬ 
persion is almost flat on a scale from 0.005 pc to 
0.1 pc (the mean velocity dispersion of molecular 
gas with /i = 2.33 is 0.204 km s _1 ), and it is close 
to the thermal velocity dispersion (0.188 km s -1 
at 10 K and fj, = 2.33). This demonstrates that 
nonthermal motions are significantly dissipated on 
scales less than 0.2 pc in the L1495-B218 filaments 
resulting in a breakdown in the cloud-scale super¬ 
sonically turbulent scaling relation. 

Observations toward molecular clouds shows 
that molecular clouds are supersonic turbulent 
(e.g. Larson 1981; Scalo 1984; Miesch & Bally 
1994; Padoan et al. 1999). In a smaller scale, ob¬ 
servations toward molecular cloud cores show that 
dense cores are typically thermally supported and 
turbulence is significantly dissipated (e.g. Caselli 
& Myers 1994; Goodman et al. 1998; Jijina et 
al. 1999; Caselli et al. 2002). The dense NH 3 
structures in this study also have the mean R p of 
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4.9 (§3.2.2). On the other hand, the R p values 
measured in C 18 0 J = 1—0 (Hacar et al. 2013) 
in the L1495-B218 filaments are close unity, and 
the mean R p value measured in 13 CO J = 1—0 
(Panopoulou et al. 2014) is 0.25. Since C 18 0 J 
= 1—0 traces less dense and larger structures than 
those seen by NH 3 (1,1) and 13 CO J = 1—0 traces 
even more extended regions than those traced by 
C 18 0 J = 1—0, the dependency of R p values on 
the effective excitation density of the tracer sug¬ 
gests that supersonic turbulence cascades from a 
cloud scale to a core scale and are considerably 
dissipated on a scale smaller than molecular fila¬ 
ments (<0.5 pc) in the Taurus. A more detailed 
analysis of the velocity dispersion and associated 
kinematic structure will be presented in Paper II. 


Mass-size relations are often pursued to investi¬ 
gate the correlation between density structures of 
sources and their size. The mass and size distri¬ 
butions of leaves and branches in the L1495-B218 
filaments are shown in Figure 14 Mass and size 


show a strong correlation of M ~ r 19 (red solid 
line for NH 3 leaves and branches) (power law in¬ 
dexes of leaves only and that of branches only are 
close to 2 as well). This is close to a global re¬ 
lation for molecular clouds where the power law 
index is 2 (Larson 1981), which indicates a nearly 
constant mean column density. Kirk et al. (2013) 
found a relation with a power law index of 2.35 in 
the Taurus from Herschel dust continuum obser¬ 
vations. Their analysis includes structures larger 
than 1 pc, whereas the largest structure in this 
study is only 0.1 pc. If we only take structures 
less than 0.1 pc from Kirk et al. (2013), the power 
law index is shallower than 2.35 and close to 2. 
Kauffmann et al. (2010) studied mass-size rela¬ 
tions resolved in several nearby molecular clouds. 
They also found that the power law index is close 
to 2 in the Taurus when size is smaller than 0.1 
pc, which agrees with our study. 


A power index of 2 suggests that the dense 
structures in the L1495-B218 filaments are grav¬ 
itationally unbound and have relatively con¬ 
stant column density while a gravitational bound 
Bonner-Ebert-like core is expected to have a power 
law index closer to 1 since p ~ r -2 for the outer 
profile of a Bonnor-Ebert-like structure. Our virial 
analysis also shows that most of NH 3 leaves in 
the L1495-B218 filaments are gravitationally un¬ 
bound. If the NH 3 leaves are formed by colliding 


flows as shown in the gravoturbulence scenario 
(e.g. Gong & Ostriker 2013), density structures 
are more likely to be shallower than p ~ r~ 2 . 

5. SUMMARY & CONCLUSIONS 

In this paper, we present extensive NH 3 (1,1) 
and (2,2) maps of the L1495-B218 filaments ex¬ 
tending over 3 degrees on the sky in the Taurus 
molecular cloud with unprecedented depth (av¬ 
erage rms is 120 mK), angular resolution (31"), 
and velocity resolution (0.038 km s -1 ). Using 
the maps, we study the physical properties of the 
L1495-B218 filaments by fitting 13109 NH 3 (1,1) 
and (2,2) spectra using an adaptive mesh refine¬ 
ment search of y 2 space. The main results are 
following. 

1. From our ammonia observations and a proto¬ 
star catalog (Rebull et al. 2010), we confirm that 
B211 and B216 are young, less evolved regions, 
while B213 and B7 are actively star-forming, older, 
more evolved regions. The young regions have 
NH 3 column density around 1 x 10 13 cm -2 and 
do not have any protostars. On the other hand, 
the more evolved, older regions have NH 3 column 
density up to 1 x 10 14 cm -2 and have multiple 
protostars. NH 3 emission is mostly confined to a 
single velocity component and the inversion levels 
are subtliermally populated. 

2. Gas kinetic temperatures in the L1495-B218 
filaments deduced from NH 3 (1,1) and (2,2) lines 
reveal that the filaments to be very cold (8 — 15 
K with median value equal to 9.5 K, 60% of spec¬ 
tra are in the range of 9.5 ± 1 K). We found that 
there is a small difference in gas kinetic tempera¬ 
ture between the more evolved regions (B7, B213, 
and B218) and less evolved regions (B10, B211, 
and B216); with the more evolved regions having 
a higher median gas kinetic temperature (by at 
most 0.5 K) than the less evolved regions. Gas ki¬ 
netic temperatures tend to decrease toward dense 
core centers at the level of a few Kelvins. This 
may ultimately affect the stability and dynamics 
of dense cores and filaments. 

3. The nonthermal velocity dispersion of NH 3 
(1,1) and (2,2) lines in the L1495-B218 filaments 
(~0.08 km s _1 ) is considerably narrower than 
those of C 18 0 J = 1—0 (~0.15 km s _1 ; Hacar et al. 
2013) and 13 CO J = 1—0 (~1 km s _1 ; Goldsmith 
et al. 2008) in the same region. This suggests 
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that nonthermal motions in the regions traced by 
NH 3 have been significantly dissipated compared 
to those found in the more diffuse portions of the 
filaments. As a result, dense regions do not display 
a Larson’s size-line width relationship. 

4. Results from the CSAR clump-finding algo¬ 
rithm found 39 NH 3 peaks (leaves) and 16 nested 
groups (branches), which is about twice the num¬ 
ber that Hacar et al. (2013) found in N 2 H+. 
This may be because NH 3 reveals lower density 
cores than does N 2 H + and our NH 3 observations 
have 67% better spatial resolution than the N 2 H + 
data of Hacar et al. (2013). The NH 3 leaves and 
branches are identified on a scale from 0.01 pc to 
0.1 pc and have masses ranging from 0.05 M 0 to 
9.5 M 0 . 

5. Masses of dense cores derived from dust con¬ 
tinuum and NH 3 show a strong correlation in the 
L1495-B218 filaments. This suggests that NH 3 
is an exceptional tracer for tracing dense cores 
in these filaments. A similar study across other 
molecular clouds will be interesting. 

6 . Most of the NH 3 leaves have observed masses 
smaller than the corresponding virial masses, 
which means they are likely to be gravitationally 
unbound structures. We found nine NH 3 leaves 
to be gravitationally bound and that 7 out of 9 
leaves are either protostellar or within branches 
associated with star formation activity, while 30 
NH 3 leaves are gravitational unbound and only 10 
out of 30 unbound leaves are either protostellar or 
within branches associated with star formation. 
We also found that 12 out of 30 gravitationally 
unbound leaves are pressure-confined. This is a 
similar conclusion that Kainulainen et al. (2011) 
found in larger scale clumps (> 0.1 pc) in nearby 
clouds, while the dense structures in this study 
are from 0.005 pc to 0.1 pc. 

7. Two sources may provide confining pressure 
to gravitationally unbound NH 3 leaves. The self¬ 
gravity of the filaments produces a gravitational 
potential well and the weight of filament materials 
outside of NH 3 leaves may exert a confining pres¬ 
sure on the surface of NH 3 leaves. The estimated 
mean pressure due to weight of filaments is a bit 
smaller than the surface pressure needed to confine 
gravitationally unbound leaves but within factor of 
two. The other source may be ram pressure due to 
large scale inflows. 12 CO and 13 CO observations 
show that there may be inflows toward the L1495- 


B218 filaments. The minimum ram pressure due 
to the inflows is about half of the surface pressure 
needed to confine gravitationally unbound leaves. 

8 . Most of the NH 3 leaves are well aligned with 
their lowest-level branches within the range of ±40 
degrees. The distribution of the apparent axis ra¬ 
tios of NH 3 leaves is closer to that of randomly ori¬ 
ented prolate spheroids than that of randomly ori¬ 
ented oblate spheroids, which indicates most of the 
NH 3 leaves are likely to be prolate. Since magnetic 
pressure would shape a dense core to be a oblate 
spheroid, magnetic fields do not seem to play an 
important role in determining the observed axis 
ratio of NH 3 leaves. 

9. The NH 3 leaves and branches each follow a 
mass radius relationship of M ~ r-i.9±o.os_ 

Finally, both core-like and filamentary struc¬ 
tures in the L1495-B218 filaments were success¬ 
fully probed using NH 3 , which suggests that NH 3 
is an excellent tracer of > 10 3 cm -3 gas for prob¬ 
ing both dense cores and surrounding filaments 
in molecular clouds. We found that most of the 
NH 3 structures agree in the LSR velocity with the 
coherent filaments found by Hacar et al. (2013). 
We also found that most of dense NH 3 structures 
(the CSAR leaves and branched) coincide with 
the dynamically critical filaments (>15 M 0 pc -1 ) 
in Hacar et al. (2013). These results suggest 
that a dense core may form as a pressure-confined 
structure, evolve to a gravitationally bound state 
within critical filaments, and then undergo star 
formation. 
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A. METHODS OF ANALYZING PHYSICAL PROPERTIES 


We derived physical quantities through finding the minimum reduced y 2 value between the observed 
NH 3 (1,1) and (2,2) spectra and the spectra model. Our method of finding the minimum y 2 value is a bit 
different from previous work in that we do not use a y 2 minimization method that utilizes a covariance 
matrix to calculate the uncertainties but instead explore a wide range of parameter space in adaptive grids 
to obtain the global minimum and determine the shape of y 2 volume as functions of the parameters. The 
common method of minimizing y 2 value utilizes a covariance matrix because the method is relatively easy 
and computationally fast, but parameters should be random variables. If a parameter has a specific non- 
random distribution, the 1 -a uncertainty of that parameter may be not symmetric and may differ from the 
standard deviation estimated by the covariance matrix. 

Exploring parameter space gives proper estimation of uncertainties, but it is usually computationally 
expensive if the parameter space is large and is searched at high resolution. To reduce computational time, 
we use the adaptive mesh refinement method (AMR method). First, we set up ranges of parameter space to 
explore [xo, lower, £ 0 ,upper], where the subscript 0 denotes zeroth layer, the subscript lower means the lower 
limit of the parameter space which we are exploring, and the subscript upper means the upper limit. The 
number of grids in each parameter can be chosen arbitrarily, but for efficiency we used fewer than 20 grids 
for each parameter. Second, we calculate the minimum y 2 & f given grid points and find its location a^min- 
Third, we set up new range of parameter space and explore as the first layer [zyiower, ®i, upper]- The range 
can be given as 


•El, lower — *E0,min 


•EO, upper *E(),lower 

2 


when lower ^ ^EO,lower 


(Al) 


and 


•El, upper — *E0,min 


«E0, upper •EO,lower 

2 


when 2Jl,upper — «E0,upper* 


(A2) 


If aq 5 lower is smaller than x 0 , lower, then aqjower = %o, lower- Likewise, if 2q, up per is larger than £ 0 ,upper, then 
ah, upper = xo, upper- Fourth, we estimate the y 2 value within the new range of parameter space with the 
same number of grids we used in the first step. At this step, the resolution is double that of the previous 
layer. Fifth, we repeat the second to the forth steps until we achieve the desired resolution. If there is the 
global minimum of y 2 value within the range of the zeroth parameter space, this code can find the global 
minimum with an uncertainty of the resolution of the final layer unless the widths of local minima and global 
minimum are much narrower than the resolution of the final layer. In our calculations, we set the finest 
velocity resolution of the model to be at least quarter of the velocity resolution of our observation, so we 
may not miss the global minimum. 

For the model of the NH 3 line profiles, we use the model described in Rosolowsky et al. (2008). 


Figure [ 15 ] shows two examples of fitting the observed NH 3 spectra and Figure [l6] shows the y 2 of the two 
examples in T k . vs. ti space. The top panel of Figure 15 is the best-fit of the simple emission model to an 
observed strong NH 3 (1,1) emission of ~32(7T mb an< i the left panel of Figure [l6| is the reduced y 2 of fitting 
in Tfc vs. t 1 space with the other three parameters fixed at the best-fit values. The bottom panel of Figure 


15 shows an example of a weaker emission spectra of NH 3 (1,1) with the peak signal-to-noise ratio (SNR) 
being ~7.5 and the right panel of Figure 16 is the reduced y 2 of fitting in T k vs. t\ space with the other 
three parameters fixed at their best-fit values. Uncertainties of physical parameters are about three times 
larger in the weak line than those in the strong line. 

We fit 13,109 spectra which satisfy the detection of NH 3 (1,1) and (2,2) emission is above 5crT mb and 
3<TT mb , respectively, where o"T mb is the baseline rms. The parameter ranges for searching for the best-fit are 
[7.5 K, 30.0 K] for T k , [2.75K,”30K] for T ex , [5 x 10" 3 , 20] for n, and [0.05 km s" 1 , 0.3 km s" 1 ] for cr v . The 
Vi sr is preliminarily estimated from the peak intensities of NH 3 (1,1) and (2,2) lines and a more accurate 
value is searched for in a range of ±0.1 km s -1 from the guessed value. Typical 1-cr uncertainties are ±0.77 
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K for Tk, ±0.85 K for T ex , ±0.26 for n, ±0.013 km s 1 for a v , and ± 0.020 km s 1 for vi sr . The overall 
uncertainty distributions of 13,109 fits to the spectra are shown in Figure [TTl 


B. MODELING APPARENT AXIS RATIO DISTRIBUTIONS OF DENSE STRUCTURES 
USING A MONTE CARLO METHOD 

We analyzed morphology of dense structures through estimating the minimum y 2 value (not reduced y 2 ) 
between the observed distribution of apparent axis ratios and distribution models. The distribution models 
are made using a Monte Carlo method. For each model, we imposed ten million spheroids and assumed all of 

the spheroids are either in oblate or prolate shapes. The apparent axis ratios of the spheroids are determined 

by the intrinsic axis ratio of the spheroid and the angle between observer’s line of sight and the equatorial 
plane of the spheroid as follow: 

p(q,6) = [1 — (1 — g 2 )cos0] 5 for oblate, (Bl) 

p(q, 9) = q[ 1 — (1 — <7 2 )sin0]5 f 01 - prolate, (B2) 

where p is the apparent axis ratio of the spheroid, q is the intrinsic axis ratio of the spheroid (short axis/long 
axis), and 9 is the angle between observer’s line of sight and the equatorial plane of the spheroid. We assumed 
that the intrinsic axis ratios of the spheroids follow a Gaussian distribution and the angles between the 
observer’s line of sight and the equatorial planes of spheroids have random values with a uniform distribution 
within [—^limiti ^limit]- The mean value and width of a Gaussian distribution and 0n m it serves as input 
parameters of the model. Using above equations with given input parameters, we modeled a distribution of 
apparent axis ratios of the spheroids. We made 400,000 models with different mean values (q) and widths 
c Tq of Gaussian distributions of the intrinsic axis ratios and - The range of (q) that is exploited in our 
models is from 0.01 to 1 with an interval of 0.01 and the range of a q is from 0.01 to 1 with interval of 0.01. 
The range of 0u m it is from 0 degree to 90 degree with an interval of 5 degree. We compared all 400,000 
distribution models to the observed distribution and found the best-fit model by minimizing the y 2 value. 
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Table 1 

Properties of NH 3 sources 
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Table 2 

Physical Properties of NH 3 sources 


ID° 

rpb 

< 

Mass from dust 

Virial Mass d 

Mass from NIT 


[K] 

[km s” 1 ] 

[M 0 ] 

[M 0 ] 

[M 0 ] 

1 

10.24^q Y 42 

n iqqU-UU 5 

U.±yy _ 0 005 

0.18 

0.45 

0.36 

1.5 

q qo 0.52 
y - yz - 0.86 

0 201 0 011 
U.ZUl_o 005 

4.56 

2.21 

3.91 

2 

9.66°q 9 08 

0 1Q7 °- 005 

u - iy ‘-0.005 

0.06 

0.19 

0.09 

3 

10.25°j 93 84 

0 204 0 014 
U.ZU4t _ 0 008 

4.38 

1.45 

1.03 

4 

8.49 9 Q 9 gg 

0 210 0 011 
U.ZIU_Q 012 

0.31 

0.34 

0.05 

5 

9.65L1 2 34 

n qc:q0-018 

u.zoz _ 0 016 

1.74 

1.93 

0.71 

6 

9.47 °_ 0 78 50 

n 9970.005 

u.zz/ _ 0 020 

0.77 

1.15 

0.50 

6.5 

9.49 9 q 9 52 

0 201 0 028 
U.ZU1_ 0 .009 

1.85 

1.61 

1.39 

7 

9.48 9 q 3 55 

n 1 Q 4 O .006 
u.iy^_ 0 .oo5 

0.88 

0.98 

0.77 

7.5 

9-29 9 q 8 77 

0 9D2 0 - 022 

u.zuz_o 009 

9.52 

3.89 

7.58 

8 

9-26 9 q 6 21 

0.196°j 9O2 O3 

0.26 

0.47 

0.37 

8.5 

8.99 9 q 2 54 

0.198°j 9O8 O5 

3.20 

1.88 

3.64 

9 

8.61°_- 30 24 

n 1 Q 7 O .007 
U.l^ / _ 0> 005 

1.02 

0.89 

1.35 

10 

9-63 2 i 2 n 

0.224°j 929 12 

0.41 

0.61 

0.25 

11 

10 . 00 °j 39 54 

0.203 9 ; 942 

0.10 

0.58 

0.12 

12 

q oqo.67 
y - zy -0.38 

f) 91 q0.026 

U.ZlO_Q 017 

0.95 

0.98 

0.98 

12.5 

9.49 °_- 94 57 

0 220 0 - 024 
U.ZZU_o 021 

1.32 

1.47 

1.24 

13 

10.90^? 6 41 

0 2 S 6 0 - 009 
u.zou_o 013 

0.04 

0.31 

0.03 

13.25 

9.68L 48 75 

n 9970.036 

U.ZZ / — 0.022 

2.79 

2.67 

2.12 

14 

9.68 9 q 9 37 

0 99fi °- 010 
U.ZZD_o 008 

0.23 

0.61 

0.13 

15 

8.84 9 o 9 50 

0.217°j 93 o 5 14 

0.21 

0.50 

0.15 

15.5 

9.50°_'f O6 

0.206°j 939 12 

2.73 

2.02 

1.88 

16 

9.66 °_ 27 30 

0 IQfiO-007 

U.l^U _ 0 003 

0.53 

0.75 

0.46 

17 

9.31°E 9 54 

0 IQO 0 - 007 
U.19O_ 0 .005 

1.62 

1.22 

1.97 

18 

9.75 9 -q 3 4o 

n 9990.006 

u.zzz _ 0 009 

0.12 

0.40 

0.09 

18.5 

9.48 9 q 4 47 

0 99D 0 006 
u.zzu_ 0 .oio 

0.99 

1.09 

0.66 

19 

9.11 °'&7 

0 9140.006 
U.Z1^_0.006 

0.07 

0.22 

0.05 

19.25 

9.45 9 ? 9 24 

0 217°- 013 
u.z± 1 _ 0 .oio 

4.06 

2.15 

2.09 

20 

9.50L- 22 29 

0 21 7 0 012 

U - Zi/ -0.007 

0.30 

0.66 

0.15 

21 

9.17^? 6 3o 

0 201 0 016 
U.ZUl_o 016 

1.92 

1.30 

0.70 

21.925 

9.80 9 q 7 82 

0.203°j 91 0 7 12 

8.97 

3.82 

9.74 

22 

9.74 9 q 8 75 


1.89 

1.39 

1.73 

22.85 

9.85°_- 82 77 

0 20 S 0 017 
u.zuo _ 0 012 

6.84 

3.55 

8.97 

23 

10.39°j 31 54 

0 210 0 005 
u.ziu_o.on 

0.95 

1.11 

1.45 

23.75 

9-91 9 ' 94 8o 

0 902 0 018 
U.ZUO_ 0 .oi 2 

4.77 

3.04 

7.16 

24 

9.08 9 q 5 16 

n 1 qi 0.003 

U.191—0.003 

0.88 

0.77 

2.14 

24.5 

9.20Lf 23 

0 1 Q 2 0 - 013 

U. 19Z_ 0 003 

1.81 

1.38 

3.63 

25 

10.26°j 48 78 

0 904 0 009 
U.ZU4t _ 0 009 

0.27 

0.53 

0.45 

26 

9.75 ^ 31 22 

0 . 202 °j 912 12 

0.24 

0.77 

0.16 

27 

10.27^ 91 og 

0.203°j 9O4 11 

0.08 

0.45 

0.05 

27.75 

10.63^ 99 31 

0.204°j 91 0 ° 14 

0.85 

1.65 

0.73 

28 

i i qoO.95 
±l.OZ _ 0 58 

0.213°j 91 o ° O5 

0.15 

0.70 

0.15 
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Table 2 —Continued 


28.5 

10.64°'f> 32 

n oncO.OlO 

U.ZUO_ 0 014 

0.77 

1.66 

0.68 

29 

10.03°y 97 

n 1 Qc;0.008 

u. iyo_ Q 007 

0.19 

0.68 

0.18 

30 

0 Q2 1 - 08 
^•^-1.18 

n 919O.008 

U.Z1Z_ 0 005 

1.08 

1.39 

0.99 

30.5 

9.77^ 4 09 

0 99f) 0 022 
U.ZZU_0.012 

3.09 

2.56 

2.03 

31 

9.55°4 9 52 

0 9S1 °- 014 
u.zoi_o.on 

0.50 

0.93 

0.29 

32 

8.95^? 6 32 

0.207°q 1 q 19 

0.26 

0.40 

0.08 

33 

9-54°j 93 58 

o^oe^on 

2.49 

1.87 

4.61 

34 

1 2 QQO.57 
lO.O^_o 91 

0.270°q 2 o 7 15 

0.07 

0.69 

0.04 

35 

9.22°? 0 n 

0.208°q° o 7 O5 

0.36 

0.60 

0.17 

36 

9-81^o 2 92 

0.205°- 31 0 2 12 

0.19 

0.57 

0.11 

37 

9.83^q° 62 

o-209™ 18 08 

1.82 

1.61 

2.97 

38 

9.55°jf 48 

n 91 9O.OO6 

u.ziz_ 0 007 

0.31 

0.62 

0.26 

39 

9.55^q 8 56 

0 iqqO.ou 

U. 19^_ 0 0O6 

2.89 

1.51 

3.22 


a IDs with decimal points are branches in dendrograms. 

k Median of the NH 3 leaves and branches. Subscripts 

and superscripts are 15.9th %tile and 84.1th %tile of the T& 
distribution within NH 3 leaves and branches. 

c Median a v of the NH 3 leaves and branches with the mean 
molecular weight /i=2.33. Subscripts and superscripts are 15.9th 
%tile and 84.1th %tile of the cr v distribution within NH 3 leaves 
and branches. 

^The NH 3 leaves and branches are assumed to be homoeoidal 
ellipsoids for the virial mass calculation. 

e The fractional abundance of NH 3 is assumed to be 1.5 X 
IQ — 9 (Typical fractional abundance ranges from 10 -9 to 10 7 , 
e.g. Jijina et al. 1999; Tafalla et al. 2004). 
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Table 3 

Geometrical Properties of NH 3 sources 


ID a 

Short axis 
[10 -2 pc] 

Long axis 
[10 -2 pc] 

Axis Ratio 

Phase Angle 1 
[Degree] 

1 

0.65 

2.27 

0.287 

38.9 

1.5 

3.79 

10.1 

0.376 

43.9 

2 

0.46 

0.8 

0.581 

44.8 

3 

2.8 

5.98 

0.469 

152.5 

4 

0.91 

1.61 

0.568 

140.1 

5 

2.75 

5.47 

0.505 

34.4 

6 

1.47 

4.21 

0.348 

151.4 

6.5 

1.67 

9.69 

0.172 

1.3 

7 

1.69 

4.85 

0.348 

18.1 

7.5 

3.37 

27.2 

0.124 

24.6 

8 

1.14 

1.99 

0.571 

20.1 

8.5 

2.21 

11.1 

0.200 

23.7 

9 

1.34 

4.57 

0.294 

16.2 

10 

3.00 

3.69 

0.813 

170.3 

11 

0.99 

1.6 

0.621 

4.3 

12 

2.09 

3.44 

0.606 

25.4 

12.5 

2.01 

5.64 

0.357 

37.8 

13 

0.64 

8.27 

0.775 

179.6 

13.25 

3.03 

11.3 

0.268 

38.9 

14 

1.16 

1.92 

0.606 

36.9 

15 

0.87 

1.82 

0.476 

38.0 

15.5 

2.86 

9.85 

0.291 

173.8 

16 

1.27 

3.67 

0.345 

143.9 

17 

4.40 

4.74 

0.935 

173.6 

18 

0.72 

1.33 

0.541 

160.2 

18.5 

2.12 

3.64 

0.581 

43.4 

19 

0.69 

0.69 

1.00 

30.0 

19.25 

5.30 

6.97 

0.763 

23.4 

20 

0.97 

2.56 

0.379 

165.5 

21 

2.67 

5.73 

0.467 

13.1 

21.925 

4.64 

20.4 

0.228 

174.4 

22 

3.81 

5.1 

0.752 

171.7 

22.85 

3.57 

21.1 

0.170 

178.9 

23 

2.11 

4.19 

0.505 

177.0 

23.75 

2.74 

19.6 

0.140 

176.8 

24 

2.25 

3.26 

0.694 

20.45 

24.5 

2.27 

7.36 

0.309 

177.5 

25 

1.13 

2.09 

0.543 

156.7 

26 

0.87 

4.31 

0.202 

4.5 

27 

0.68 

2.09 

0.328 

41.2 

27.75 

2.6 

7.41 

0.352 

20.7 

28 

1.82 

2.32 

0.787 

37.8 
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Table 3 —Continued 


28.5 

2.16 

8.22 

0.263 

16.4 

29 

1.19 

3.32 

0.361 

4.2 

30 

1.94 

5.96 

0.326 

142.3 

30.5 

2.47 

12 

0.207 

147.9 

31 

1.52 

2.88 

0.529 

137.7 

32 

1.19 

1.38 

0.870 

138.8 

33 

3.58 

7.46 

0.481 

7.0 

34 

0.84 

1.57 

0.538 

159.2 

35 

1.29 

2.26 

0.571 

27.3 

36 

1.01 

2.41 

0.418 

28.3 

37 

3.4 

5.98 

0.571 

2.4 

38 

1.28 

2.23 

0.571 

17.5 

39 

4.45 

5.79 

0.769 

146.8 


a IDs with decimal points are branches in 
dendrograms. 

^The phase angle of NH3 leaves with re¬ 
spect to the Galactic longitude. For the 
phase angle with respect to the equatorial 
plane. 62 . 6 ° should be added. 
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Fig. 1.— Top: 500 /rm dust continuum emission (color) seen by the SPIRE instrument of the Herschel 
Space Observatory. Bottom : map of integrated intensity of NH 3 (1,1). 
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Fig. 2.— rms of the ammonia T m b map at a spectral resolution of 6.1 kHz (corresponding to 0.077 km s 1 ), 
obtained by smoothing two spectral channels. 
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Fig. 3.— Maps of integrated NH 3 (1,1) emission and positions of known protostars. X-axis is galactic 
longitude and Y-axis is galactic latitude. The names of previously studied cores are indicated in the image. 
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Fig. 4.— Map of the gas kinetic temperature in L1495-B218 filaments. X-axis is galactic longitude and Y- 
axis is galactic latitude. Locations of protostars are marked with various symbols depending on the spectral 
classification of protostars. Red arrows indicate deg®e cores having gas kinetic temperatures decreasing 
toward core centers. Regions with multiple velocity components in B211 are suppressed. 












Fig. 5. Histograms of physical properties for each region. Vertical lines denote the median values. High 
fractions at the very left edge of the first panel is due to low signal-to-noise level at the edge of map. Those 
are excluded in estimating median values. B211 is also excluded due to the presence of multiple velocity 
components. The dotted vertical line in the fourth panel denote R p = 1. 
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Fig. 6.— NH 3 (1,1) lines (black histograms), the best-fit single component NH 3 (1,1) model (dashed red 
lines), C 18 0 1—0 (green histograms; Hacar et al. 2013), and 13 CO 1—0 (Goldsmith et al. 2008; Narayanan et 
al. 2008) at the NH 3 peak positions of NH 3 leaves identified by CSAR. Black vertical lines denote locations 
and relative strength of NH 3 (1,1) hyperfine lines. Red-green arrows denote the LSR velocities of NH 3 leaves. 
13 CO and C 18 0 have multiple velocity components in some of NH 3 leaves but NH 3 seems to be confined to 
one velocity component and is well fitted with a NH 3 line model with a single velocity component. 
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Fig. 7. — Continued from Fig. [6] 
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Fig. 8.— Velocity dispersions of nonthermal components. X-axis is the galactic longitude and Y-axis 
is the galactic latitude. Locations of protostars are marked with various symbols depending on spectral 
classification of protostars. Regions with multiple vejplcity components in B211 are suppressed. 











Fig. 9.— Location of NH 3 leaves. Solid and dashed orange contours are NH 3 leaves and branches, respec¬ 
tively, identified by CSAR. 
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Fig. 10. Dendrogram of NH 3 sources. Red solid lines denote active leaves (leaves that are embedding 
a Classl/Flat spectra protostar within their lowest-level branches or closely associated with a Classl/Flat 
spectra protostar within three times their leaf radiuses) and dashed red lines denote active leaves with an 
embedded protostar. Black circles indicate mass of NH 3 sources estimated from 500 fj,m dust continuum. 
Red circles denote the virial mass. Purple dashed vertical lines indicates different subregions. 
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Fig. 11.— Location of dust and NH 3 leaves. Blue contours are boundaries of dust leaves and green-black 
crosses are the peak positions of NH 3 leaves. 


37 












Moss from Dust [M Sun ] Virial Moss [M 





Fig. 12.— Relations between masses (the first, second, and third panels). The red line denotes x=y. Ratio 
of thermal support to nonthermal support (the fourth panel). The red line denotes R p = 1. 
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Fig. 13.— Distribution of properties of the NH 3 leaves. Red histograms denote active leaves, and blue 
histograms are inactive leaves. The first panel is the ratio of the virial mass to the observed mass. The red 
and blue vertical lines are the mean ratios of active and inactive leaves, respectively. The second panel is the 
ratio of the work done by pressure to the gravitational energy at the surface of NH 3 leaves. The red and blue 
vertical lines are the mean ratios of active and inactive leaves, respectively. The black dashed vertical line 
is the energy ratio of the critical Bonnor-Ebert sphere. The third panel is the angle difference of the NH 3 
leaves to their branches. The fourth panel is the axis ratio of the NH 3 leaves deduced from the principal 
axis analysis. The solid and dashed lines are the best-fit probability density of apparent axis ratio of prolate 
and oblate spheroids with a uniform orientation angle distribution within [-80°, 80°], respectively. 
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Fig. 14.— Relations between physical properties of NH 3 sources. The left panel shows the relation of the 
velocity dispersion (/z=2.33) and the size of the NH 3 leaves. The red dotted line denotes the thermal velocity 
dispersion when T^ = 10 and /i = 2.33. The right panel is the mass-size relation of the NH 3 sources. Black 
filled circles are the NH 3 leaves and blue boxes are the NH 3 branches. The red solid line is our best fit of 
M~ r 19 to the NH 3 sources. The purple dashed line is M~ r 1 - 5 from Kauffman et al. (2010) and the orange 
dashed-dotted line is M~ r 2,35 from Kirk et al. (2013). 
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Fig. 15.— Examples of spectra model fitting (red dashed lines) to observed lines (black solid lines). The top 
panel shows NH 3 (1,1) lines having a peak signal to noise ratio (SNR) equal to 32, and the bottom panel 
shows NH 3 (1,1) lines having a peak signal to noise ratio equal to 7.5. 
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Fig. 16.— Examples of reduced the % 2 space in vs. t\. The left panel shows the x 2 space of fitting the 
line with a peak SNR of 32, and the right panel shows x 2 space of fitting the line with a peak SNR of 7.5. 
The white contours are l-er uncertainty limits. 
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Fig. 17. Distributions of l-cr uncertainties in the specti 
values. 
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i. Vertical lines represent the median 






